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Abstract 

The thermodynamic second law in the evaporating black hole space-time is ex- 
amined in the context of two-dimensional dilaton black hole. The dynamical evo- 
lution of entropy is investigated by using the analytical perturbation method and 
the numerical method. It is shown that the thermodynamic second law holds in the 
vicinity of infalling shell and in the region far from the infalling shell in the case 
of semi-classical CGHS model. The analysis of the general models is also presented 
and its implication is discussed. 



1 Introduction 



At present, the thorough understanding of quantum gravity has not been achieved yet. To 
reveal the nature of quantum gravity, it is important to figure out the quantum aspects of 
black hole. The existence of the non-trivial causal space-time structure such as black holes 
gives a profound insight into the character of gravity. Curiously, the geometry of black hole 
is intimately related to thermodynamics The discovery of Hawking radiation stressed 
the role of the black hole thermodynamics]^]. It also gives the picture of an evaporating 
black holes. Recently, the problem of information loss for the evaporating black hole has 
been discussed widely. Originally, the thermodynamic laws of black hole were proved in 
a stationary situation. Hawking's calculation was performed on a fixed stationary black 
hole. However, Hawking radiation naturally leads us to the evaporating black hole picture. 
Hence, it is legitimate to ask how the thermodynamic laws are modified. 

In this paper, we shall investigate the thermodynamic second law in the context of 
two-dimensional dilaton gravity. The two-dimensional CGHS model ^ provides an inter- 
esting toy-model for the study of black hole thermodynamics. For this model, it is easy to 
give the classical general solution. There exists a black hole solution. Moreover, the ther- 
modynamic laws hold in this model]^ and the standard calculation yields the Hawking 
radiation with a certain temperature. It was shown that the effective action in which the 
backreaction of the Hawking radiations is included gives the evaporating black hole pic- 
ture. On the other hand, Wald presented a new way to look at the entropy of black hole as 
a Noether charge associated with diffeomorphism invariance of a theory]^ [§. In particu- 
lar, the first law of black hole mechanics is a general relation between "boundary terms" at 
event horizon and spatial infinity in a stationary asymptotically fiat space-time. Iyer and 
Wald derives an inductive algorithm for general diffeomorphism invariant theories 0, and 
the expression of the entropy is given by 



in the context of two-dimensional dilaton gravity. Here L{gab, R) is the Lagrangian of 
the theory and Qab is the metric of the spacetime and R is the scalar curvature. Myers 
has proposed to use the above entropy expression in dynamical situations and has shown 
that the entropy satisfies the thermodynamic second law in the classical and semi-classical 
RST model ^ . However, RST model is too special to conclude that the thermodynamic 
second law holds in general. The question we would like to address is whether the Myers's 
result can be extended to unsoluble semi-classical CGHS model and the more general 
models. 
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2 Entropy of semi-classical CGHS model 
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2.1 Classical model and its solution 

To make this paper self-contained, we start from classical CGHS model. The classical 
CGHS action is 

1=^1 d'xV^g{e-"f[R + 4(y<pf + - - ElV/.)^ (2) 

where (j) is dilaton field, R is 2d Ricci scalar, A is a positive constant, V is the covariant 
derivative associated with 2d metric g^^i^, and fi are N massless scalar fields. In conformal 
gauge, ^ 

9+- = "2^^''' = 9— = ^, = t± r, 
equations of motion for p, 0, fi are 

e-''^(2 9+9_(/)-49+(/)9_(/)- A'e''') = 0, (3) 

-4 9+9_(/) + 4 9+09_(/) + 2 9+a_p + A^e^'' = 0, (4) 

d+d-fi — (for each i), (5) 

with the constraints 

1 ^ 

e-2'^(4a±pa±0 - 2dl<f>) + -Y^d^fid^fi = 0. (6) 

i=l 

As we can take (p = p gauge for this model, exact solution of this model can be obtained 
easily. Vacuum solution of this model in = p gauge is 

e-^P = e-^'^ = M - X'^x+x-. (7) 

It can be easily shown that the Penrose diagram of this space time is the same as that of 
Schwarzschild's spacetime. This fact shows that CGHS model has a black hole solution. 
When M — 0, the solution becomes 

^-24> ^ g-2p ^ _x''x^x-, (8) 
and we have the fiat space time 

= -da+da-, (10) 

where we have introduced a coordinate system by Ax^ = ie^*^*. The solution is called 
linear dilaton vacuum(LDV). 

Let us consider the matter falling along the null line, x~^ = Xq, then the energy- 
momentum tensor of the scalar fields becomes 

1 M 

-d+fid+fi = -y-^K^'^ - ^0 )' 

\d.f,d.f\ = l^+f^^.fi = 0. 
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The solution for this is 

^-2<t> ^ g-2p ^ - Xt)e{x+ - X+) - A'x+X", (11) 

where M can be shown to be the ADM mass of black hole. 

2.2 Evolution of entropy 
2.2.1 Entropy formulation 

Now we turn to the quantum theory. Though we cannot perform the full quantization 
of gravity, we can calculate matter contribution and obtain effective action. In covariant 
non-local form, quantum correction which is called Liouville action is 

h = ~^j d'xV^R^^R (12) 

where k = N/12, N is number of scalar fields. 

Using Wald's scheme mentioned in the Introduction, the contribution to entropy from 
the classical action is 

So = 2e-^t (13) 

where the right hand side is evaluated on the horizon. As for the non-local action, careful 
argument is necessary in the Noether charge analysis. Myers showed that such a non-local 
term actually can be handled The result is 



Si = -Kj d'y^-g{y) G{xh, y)R{y) (14) 

where G{x, y) is the Green function, that is, 

VlG{x,y)=5\x,y) (15) 

We have to specify the boundary condition for the Green function G(x, y) in order to 
make the expression (0) have a definite value. Combining (|13|) and (0), we obtain 

S = 2e-^^ + 2KP + K log{-\'^x+x-), (16) 

in Kruskal coordinate, and the last extra term came from boundary condition of (p!5|). 

The entropy (0) has a local geometrical expression and is to be evaluated at a certain 
part in spacetime. In stationary spacetime, it is the killing horizon. In dynamical situ- 
ation, on the other hand, it seems to be reasonable to consider the entropy on apparent 
horizon since it is physically meaningful and a locally-determined structure formed in the 
djTiamical spacetime. Following this proposal, we address the question of the validity of 
the entropy increase along the apparent horizon in this paper. 

Surely we have to know the solutions of and p in order to evaluate above entropy 
formula. Unfortunately, however, we cannot solve the equations of motion in the semi- 
classical theory. Thus we analyze the behavior of entropy in the vicinity of the infalling 
shell at which apparent horizon starts. After that, some results of numerical calcula- 
tions will be presented because the above analytic calculation cannot give the behavior of 
entropy far from the infalling shell. 
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Figure 1: This figure shows schematic explanation of situation considered here. 



2.2.2 Analytic approach 

We would like to evaluate the entropy along the apparent horizon. Fortunately, if we 
restrict our consideration to the vicinity of the infalling shell (Figure 1.), one can calculate 
the variation rate of the entropy along the apparent horizon as 



dS 

dx^ 



d+S + 



dx 



dx~^ 



2Kd+p + 



dx+ 



(17) 
(18) 



where we used the condition for the apparent horizon (9+0 = 0), and all of the quantities 
are evaluated at the shell|jTO|. From the equations of motion|jTl|, one can show that 



dx~ 



dx^ 



So, we must have (9+p, 9^0 and 9+9_0 at the intersection point between the apparent 
horizon and the shell line. 

For the semi-classical action (H) (|12|) , the equations of motion are 



-20 



(4^4 



4(9+ 9_0 + 2(9+ (9_p + X^e^P) = 0, 



(19) 



e-2'^(49+09_0 - 29+9_0 + A^e^O + Kd+d^p = 0. (20) 

On the shell line (x^ = Xq), we define S = 9+0, S = 9+p and w = e~'^'^ = e~^^ = 
— A^x^x", which is the LDV solution. Boundary conditions for and p is that they 
asymptotically coincide with the classical solution because quantum correction would be 
negligible in the asymptotic region. 

With these variables, (p!9[) and (^) can be written as 

{2w - k) du,E - 2wd^T. = 0, (21) 
^a^S - 2wd^J: - 2S - ^ = 0, (22) 

Xq 
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where we have used 9„0 = X^Xq /2w^ 9_ = —X^x^dw By eUminating dj^, we have 

1 1 

Co 



AL' + -A'S = -^A', (23) 



where ' denotes dw, and A, A' are defined as A = 4w(ty — k), A' = 5^y4 = — An. This 
differential equation can be solved easily. With the boundary condition, we have 

„ 1 M 1 



2xq Xx^ 

Because the apparent horizon is defined by d+cf) = S = 0, we can find the point from 
which the apparent horizon starts. Then we have 



A 







-.^i = 4f^+(-) (25) 
or in terms of w. 



We can find other solutions which satisfy (|25D. They are, however, not interesting because 
they are behind singularity. More details will be mentioned later. 

With the above solution of E, S is obtained by integrating equation 



Using these results, we have 



2M , /2M 



Next, we calculate the other part of dS/dx^. The gradient of the apparent horizon 
contains d+d^ip and 9^0, the former can be estimated from the equations of motion 
directly. Eliminating d^d^p from (|T^) and (^Dj) one has the value of d^d^cf) on the apparent 
horizon as 

^+^-<^U.H. = Y^A^. (29) 

Now, we define a new variable JF = d'^cj). Differentiating (|19D and ( ^Of ) and eliminating 
d\p, we have the differential equation for 

1 A' 
AT' + -A^T = A'd^iwi:'^) - —^{E - S) + 4kwS(9^(S - S). (30) 



Solution of this differential equation with boundary condition is 



2x, 



1 M 
— + 





M A! M2 
+ 



M2 A' 2M 2M2 1 
_ + 



2k\^xI a \xl 



3kM A' 



+ 



+ 



M 



(31) 



Substituting and A'^ into above result, one finds the value of T at the apparent horizon, 
and finally we have 



dS 




'4M2 



4M2 , 
+ Tir-rlog 



X^x 







2M 



24 



K 



(32) 



It is easily seen that the last log term is always positive. It can be proved that the sum 
of the rest terms are always positive. Therefore we could show that the entropy does not 
decrease along the apparent horizon at least in the vicinity of the infalling shell. 



2.2.3 Numerical approach 

To know the behavior of entropy along the apparent horizon far from the shell, we must 
perform numerical calculation for some values of the parameters k and M. Although we 
cannot cover all values of the parameters in numerical calculation, it is necessary to know 
information which cannot be obtained from the above analytic approach. 

Numerical calculation was performed on coordinate, and we fix the parameters, 
A = 1, (Tg" = {xq = l)^^. Equations for (p and p are second order partial differential 
equations. We have used the numerical schemes invented in |]T3 . The range of is from 



the infalling shell to the merging point of the apparent horizon and singularity. One of 
the typical numerical results is shown in Figure 2. The result is for n = M = 1. Solid 
curve is numerically calculated entropy, and dashed strait line is analytically calculated 
one. In the latter one, the value of entropy on the infalling shell is evaluated by using 
the LDV solution, and gradient of entropy is calculated with analytic scheme presented 
in above subsection. This result shows that the entropy grows still even in the region far 
from the infalling shell. The range in which analytic approximation is good is not so wide 
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Figure 2: This figure shows the entropy of semi-classical CGHS model for parameters 
A = 1, = (xq =1), K = M = 1. Solid curve is numerically calculated entropy, and 
dashed strait line is analytically calculated one. 



in this case, it shows the higher order effects make the growing rate more rapidly. We 
also calculated the evolution of the entropy for other parameters. All of the results show 
the same feature. 



3 Entropy of generic 2 dimensional model 



3.1 Generalized action and entropy 

The semi-classical CGHS action is given by adding the non-local counter term Ji. Other 
covariant local counter terms can be added to the action because of ambiguity in defining 
measures used in the functional integral, and combinations of these counter terms make 



different models [1131. Possible local counter terms which can be added to the action are 



^ / d^x^{a^R + P{V(Pf}, 

OTT J 



(33) 



K 

8^ 



/ rf\v^f;{a„0"i? + 6„0"-i(V0)2}. 

n=2 



(34) 



Because we will consider black hole solutions in these models, which are asymptotically 
flat, it is reasonable to require that models include the LDV solution. This requirement 



leads ar, 



0, so the action h will not be considered below. If 1 



(3 Will noL ue consiuereu ueiow. ii±— -g— 'j = 0, the 
model is exactly soluble because of the symmetry which leads d+d-{p — 0) = just like 



the case of classical CGHS model [ITl 



Equations of motion derived from the generalized action I = Iq + Ii + I2 are 



e-2*(49+0 (9_0 - 49+a_0 + 2d+d-p + A'e"^ + '-{ad+d.p + (3d+d^ 



(35) 



e-^'^(4a+09_0 - 29+9_0 + A^e^O + K{d+d^p - 



(36) 
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and the constraints are 



1 ^ 
^ i=l 

+«:{(5±p)' - dip + t± - ^(2a±0 9±p - dl<P) - ^{d^<Pf} = 0. (37) 
The determinant of the coefficient matrix of the kinetic term is 



In the case 1 — a/2 — [3/4 < 0, there is no degeneracy because the above formula always 
takes positive value. It means that the solution has no singularity. As we would like 
to consider the entropy of the evaporating black holes, this case is excluded from our 
consideration. By substituting the LDV solution (H) into the constraint (|37D, the boundary 
terms are determined as 

If P = 0, this corresponds to the boundary conditions adopted in the CGHS model and 
RST model. We then have 

< >= -f — , (40) 
which are evaluated from the Liouville part of the constraints (|37D. This means that the 



LDV spacetime is ffiled with radiation ||T5[. We consider only the case /5 < as physical 
system. 

3.2 Analytic approach for generic models 

With these general action and equations, we can also investigate the behavior of entropy 
in generic model. Extra term of entropy which did not appeared in CGHS model is 

S2 = -^ct>, (41) 

that comes from the local counter term I2. Thus the resulting formula for the entropy is 
given by 

S = 2e-2^ + 2Kp~^(f) + K \og{-X^x+x-). (42) 

Calculation of the entropy evolution near the in-falling infalling shell is almost similar 
to the case of semi-classical CGHS model. Equations of S and S are 

|2w; + K (^^ - 1^^ d^E - |2u; - ^(a + /3)| = 0, (43) 

Kdu,^ - (2w + ^) (9^S - 2S - 4r = 0. (44) 
V 4 / Xq 
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The solutions are 



1 M 1 

"24 ^ Xxt Vb' 



(45) 



and 



— (B' + Ak)J: + ^—B'--^, (46) 



where B and 5' are defined to be 



B' = 4|2i« + k(^-i)|, (47) 

and these can be reduced into A and A' ii a = f3 = 0. 

Solving E = 0, one obtains the point from which apparent horizon starts. For 1 — 
a/2 — (3/ A > 0, there are two solutions. We take the largest one as a solution, wq, because 
the other one is hidden behind of the singularity. By solving S = 0, 



fin 



2M 



With the above results, contribution of d+S at apparent horizon can be evaluated. 

Contribution of d^S can be calculated in the same way to CGHS model. Differential 
equation for has also the same form as (pO|), 

5jr' + Ib'T = B'd^{wT?) - -^(S - S) + 4/€u;S9^(S - S). (49) 

We must consider two cases separately, one is 1 — a/2 — /3/4 = 0, another case is 1 — 
a/2 — /3/4 > 0. Semi-classical CGHS model belongs to the latter case. 



3.2.1 Exactly soluble case 



In the case of 1 
Bose,et.a/ 



a/2 — P/4: = 0, which is exactly soluble like RST model or model of 
= p gauge, so that the latter two terms of (H) are zero. 
Thus the solution and the value at the apparent horizon can be found easily. Then we 
have the time evolution of entropy in the exactly soluble case as 



T5|] , we can take 



dS 



K,a 



2M 
X 



2M 
A 



4f-l)' 



(50) 
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This takes positive value when model-fixing parameter satisfies 



or equivalent ly, 



;2M 
0<„<4| — + 



In the case of RST model (a = 2, /9 = 0)[p|, the increasing rate of entropy is 



dS 



K 



2AI 



— > 



and in the case of Bose,et.a/ model (a = 4, f3 = —4), 



dS 
dx~^ 



K 

= — >0. 

AH. Xq 



(51) 



(52) 



(53) 



(54) 



3.2.2 Unsoluble case 

In the other case, 1 — | — | > 0, the solution of (|^) cannot be reduced into one of exactly 
soluble case, so we must distinguish this case from the previous one. With the solution of 



and contribution of S+S", we have 



dS 



dx^ 



n— ) 




+ 



5^-4«(f -1) 
32M3 1 2^2 



/t^A^Xq Xq 

2 f 1 - f - ^ 







1 _ " _ ^ 
2 4 



4M2 
kX^Xq 



X^X^n 



47 + 1-1 5o + Wl-f-| 
-log 



1 _ £ _ ^ 

-^2 4 



B'o 



AkJI - f - f 



(55) 



(56) 



The sign of this value depends on a and j3. For example, in the case a = 12, /? = —20.1 
at A = 1, Xfj" = 1, which has slightly different parameters from that of the exactly soluble 
case, dS/dx~^ takes the negative sign, leads also negative dS/dx~^ near k = M. Therefore 
there exist some models which have decreasing entropy with Wald's scheme in 2 dimension. 



4 Conclusion 

In this paper, we have investigated the dynamical evolution of the entropy proposed by 
Wald along the apparent horizon in the context of the 2-dimensional dilaton gravita- 
tional theory. In particular, we have considered the quantum backreaction and then the 
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evaporating spacetimes. In the case of the semi-classical CGHS model, we have shown 
analytically that the thermodynamic second law holds along the apparent horizon in the 
vicinity of the infalling shell. Furthermore, the numerical calculation is also presented, 
which suggest that the thermodynamic second law holds even in the region far from the 
infalling shell. Myers's result on the RST model and analysis on the semi-classical CGHS 
model are very impressive, however, they are yet special models. Hence, we have also 
studied the general models which seem physically natural in our judge. Analytic calcu- 
lation shows that, even in the vicinity of the infalling shell, the thermodynamic law does 
not hold in certain models. There are two possible choices here. One is to use this fact 
for restrictions of models. Another is to modify the entropy formula. In the latter case, 
the statistical interpretation of the entropy would be necessary. However, this interesting 
subject is beyond the scope of this paper. Further implications of our analysis should be 
considered in the future. 
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